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Abstract 

The phase transition between the valence-bond-solid (VBS) and nematic 
phases, the so-called deconfined criticality, was investigated for the quan- 
tum S = 1-spin model on the spatially anisotropic triangular lattice with the 
biquadratic interaction by means of the numerical diagonalization method. 
We calculated both VBS- and nematic-order parameters, aiming to clarify 
the nature of this transition from complementary viewpoints. Simulating 
the clusters with N < 20 spins, we estimate the correlation-length critical 
exponent as v = 0.95(14). We also calculated Fisher's exponent (anomalous 
dimension) for each order parameter. 
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1. Introduction 

According to the deconfined-criticality scenario [l|, 0, Eij], in two dimen- 
sions, the phase transition separating the valence-bond- solid (VBS) and anti- 
ferromagnetic phases is continuous. (Naively j2[ , such a transition should be 
discontinuous, because the adjacent phases possess distinctive order param- 
eters such as the VBS-coverage pattern and the sublattice magnetization, 
respectively) Extensive computer simulations have been made to support 
this claim |H, However, still, one cannot exclude a possibility of a weak- 
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first-order transition accompanied with an appreciable latent-heat release 

0, H B B E3, Q H 0, Q. 

The magnetic frustration is a good clue to the realization of the VBS 
phase 14, 15, 16j]. Alternatively, one is able to stabilize the VBS state with 
the spatial anisotropy and the biquadratic interaction [l7|, [l8| . In our pre- 
ceding paper 19l . we investigated the quantum S = 1-spin model on the 



spatially anisotropic triangular lattice with the biquadratic interaction, Eq. 
(QQ), and analyzed the singularity of the VBS-nematic phase transition; a 
schematic phase diagram is presented in Fig. (TJ where the parameters J 
and j control the spatial anisotropy and the biquadratic interaction, respec- 
tively. Scrutinizing the scaling behavior of the excitation gap, we obtained 
the correlation-length critical exponent v = 0.92(10) [19J. In this paper, we 
calculate both VBS- and nematic-order parameters, and survey the criticality 
from complementary viewpoints. 

To be specific, we present the Hamiltonian for the S = 1-spin model on 
the spatially anisotropic triangular lattice with the biquadratic interaction 



(i) 



Here, the quantum S = 1 spins {Sj} are placed at each triangular-lattice 
point i; see Fig. |5](a). The summation YjUj) runs over an possible 

nearest-neighbor (skew-diagonal) pairs. The parameter J(> 0) (J') denotes 
the corresponding coupling constant. Hereafter, we consider J' as the unit 
of energy (J' = 1). Along the J bond, both quadratic and biquadratic 
interactions exist, and the parameter j(> 0) controls a strength of the former 
component. The J'-bond interaction is purely biquadratic. The interaction 
J interpolates the one- dimensional (J = 0) and square-lattice (J — > oo) 
structures; correspondingly, there appear the VBS 20j and spin-nematic [2l( 
phases (Fig. [T]). In order to take into account such a geometrical character, 
we implement the screw-boundary condition [Fig. [2] (b)] through resorting 
to Novotny's method [ij (Section EZQ. 

The critical indices for the deconfined criticality have been investigated 
thoroughly. We overview a number of related studies. First, the 5 = 1- 
spin model on the spatially anisotropic square lattice with the biquadratic 
interaction [17| was simulated by means of the quantum Monte Carlo method; 
here, the spatial-anisotropy axis is set to be parallel to the primitive vector 
of the unit cell. At the VBS-nematic phase boundary, the authors estimated 
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the (reciprocal) correlation-length critical exponent as yv[— ( z/l/ ) -1 ] — 2.8(2) 
and Dq[= (i^) -1 ] = 2.8(2) for the VBS- and nematic-order parameters, 
respectively. These indices coincide, as anticipated; namely, a relation v = 
v v — v q holds. Additionally, they obtained Fisher's exponent for the nematic 
order t\q = 1.25(20) — z\ postulating the dynamical critical exponent z — 1 
(see below), one arrives at t/q = 0.25(20). (Its VBS-order counterpart rjy 
was not given.) In our model, Eq. ([TJ, we impose the spatial anisotropy 
along the skew-diagonal direction [Fig. 12(a)], aiming to make an alternative 
approach to this problem. Second, the internal symmetry of each spin is 
extended from an ordinary SU(2) to SU(iV) with an arbitrary integral number 



N [241 . |25| . l26l |27| or even a continuously variable parameter N [28] ; note that 
the S = 1-spin model with the (finely-tuned) biquadratic interaction can 
possess the SU(3) symmetry, and this problem is quite relevant to ours. The 
latter model |28j demonstrates a clear evidence of a deconfined criticality with 
the critical indices, 0.75 < v < 1, z — 1, and r\ = 0.63(4) at N c = 4.57(5); 
here, the index r\ denotes the ordinary Fisher's exponent for the constituent 
SU(iV) moment. It would be noteworthy that Fisher's exponent acquires 
a considerable enhancement. Third, pioneering considerations jH, [ij on the 
deconfined criticality were set forward for the S = 1/2 square-lattice model 
with the biquadratic (plaquette-four-spin) interaction. Sets of critical indices, 
v = 0.78(3), v = Vv = 0.26(3), and z = 1 fl; v = 0.68(4), rj = 0.35(3), and 



z = l|5j, have been reported. A closely related model, the so-called Q v model 



29|, yields almost similar results, v = 0.78(3), (3 = 0.27(2), and f3 v = 0.68(3); 
here, the indices, /3 and /3y, are the magnetization critical exponents for the 
constituent moment and the VBS-order parameter, respectively. Last, we 
recollect a number of field-theoretical considerations. A formalism for the 



S = 1-spin model is provided in Ref. [30(, where a microscopic origin of 
an enhancement of Fisher's exponent is argued. An extensive simulation on 
the hedgehog-suppressed 0(3) model [3l| revealed a clear evidence for the 
novel critical indices, v = 1.0(2) and r] = 0.6(1). As a reference, we quote 
the critical exponents for the three-dimensional Heisenberg universality class 



32|, v = 0.7112(5) and 77 = 0.0375(5). 

The rest of this paper is organized as follows. In Section HI we present 
the simulation results. Technical details are explained as well. In Section |3l 
we address the summary and discussions. 
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2. Numerical results 



In this section, we present the simulation results. Before commencing 
detailed analyses of criticality, we explain the simulation algorithm in a self- 
contained manner. Throughout this section, we fix the parameter j to an 
intermediate value, j/J' = 0.5, where the finite-size-scaling behavior im- 



proves [19(|; see the phase diagram, Fig. [TJ The number of spins (system 
size) ranges in iV = 10, 12, ... , 20. The linear dimension L of the cluster is 
given by 

L — y/N, (2) 

because N spins form a rectangular cluster. 

The Hamiltonian ([I]) possesses a number of symmetries, with which one is 
able to reduce (block-diagonalize) the size of the Hamiltonian matrix. Here, 
aiming to eliminate the Hilbert-space dimensionality, we look into the sub- 
space with the total longitudinal-spin moment J2iLi &i = 0> even parity, 
and the wave number k = with respect to the internal-spin-rotation-, spin- 
inversion- (S£ —> — Sf), and lattice-translation-symmetry groups, respec- 
tively. The ground state belongs to this subspace. The size of the reduced 
(block-diagonal) Hamiltonian is 9436203 for iV = 20. 

2.1. Simulation method 

As mentioned in the Introduction, we impose the screw-boundary condi- 
tion on a finite cluster with N spins; see Fig. [2](b). Basically, the spins, {Si} 
(i < N), constitute a one-dimensional (d = 1) structure, and the dimension- 
ality is lifted_to d = 2 by the bridges over the long-range pairs. According 
to Novotny [22] , the long-range interactions are introduced systematically by 
the use of the translation operator P; see Eq. (j^J), for instance. The operator 
P satisfies the formula 

P\S 1 , S 2 , ■ ■ ■ , Sn) = \Sn, Si, ... , Sjv-i)- (3) 

Here, the base diagonalizes each of {Sf}] namely, the relation S^\{Si}) = 

Sk\{Si}) holds. Novotny's method was adapted to the quantum S = 1 XY 



model in d = 2 dimensions 23 . Our simulation scheme is based on this 



formalism. In the following, we present the modifications explicitly for the 



sake of self-consistency. The XY interaction Hxy, Eq. (4) of Ref. 23], has 
to be replaced with the Heisenberg interaction 

N 

H xxx (v) = J2(P V S?P~ V S? + P v SfP" v S^ + P V S*P- V S*). (4) 

i=l 
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Additionally, we introduce the biquadratic interaction 

N 5 

2 x " ' 2 



N 5 

= -.few + 5 X) E P v Q?P~ v Qf- (5) 



Here, we utilized an equality 



(S, ■ S,) 2 = -S, ■ S s /2 + £ Q?Q? /2 + 4/3, (6) 



Q=l 



with the quadrapole moments, Qj = (Sf) 2 - (Sf) 2 , Qf = [2{S'tf - {Sff - 
(^i) 2 ]/ v^3, Qi = SfSi + SfSf, Qf = SfS* + S?Sf, and = SfSf + SfSf. 
Based on these expressions, we replace Eq. (3) of Ref. [23j with 



H = -J[jH xxx (VN)+jH xxx (VN-l)+H 4 (VN)+H A (VN-l)]-J'H 4 (l). 

(7) 

We diagonalize this matrix for N < 20 spins. The above formulae complete 
the formal basis of our simulation scheme. However, in order to evaluate the 
above Hamiltonian-matrix elements efficiently, one may refer to a number of 



techniques addressed in Refs. [22l . |23 



2.2. Critical behavior of the VBS-order parameter 

In this section, we investigate the critical behavior of the VBS-order pa- 
rameter; in the VBS phase, along the J'-bond direction, the staggered-dimer 
order develops (2o| . 

In Fig. |3l we present the Binder parameter 



U V = 1 - — (8) 



for various J(/ J') and N = 10,12, ... , 20. Here, the quadrature of the VBS- 
order parameter is given by 

<=^EW^ + A% (9) 

i=2 

and the symbol (■ • ■) denotes the ground-state average. The interaction pa- 
rameter j is set to an intermediate value j = 0.5, as mentioned above. As 
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the system size increases, the Binder parameter increases (decreases) in the 
long- (short-) range-order phase. Hence, in Fig. EJ we see that the VBS order 
develops in the small-J regime; low- dimensionality promotes the formation 
of the VBS state. The intersection point of the curves indicates a location of 
the critical point. 

In order to extrapolate the critical (intersection) point to the thermody- 
namic limit, in Fig. HI we plot the approximate critical point J^(Ni,N 2 ) 
(pluses) for [2/(Ni + N 2 )] 2 with 10 < Nt < N 2 < 20; the parameters are the 
same as those of Fig. |3j Here, the approximate transition point, J^(Ni, N 2 ), 
denotes a scale-invariant point with respect to a pair of system sizes (Ni, N 2 ). 
Namely, the relation 

Ua{Ni)\j = j ? ( NuN2 ) = U a (N 2 )\j=J ? (N 1 ,N 2 ), (10) 

with a = V holds. The least-squares fit to the data of Fig. H] yields an 
estimate = 0.236(18) in the thermodynamic limit, N — > oo. A remark 
is in order. The wavy character in Fig. H] is an artifact due to the screw- 
boundary condition. Actually, there appears a bump (drop) at N « 4 2 (« 3 2 , 
5 2 ); in other words, an undulation comes out, depending on the condition 
whether the linear dimension of the cluster, y/~N, is close to an even number 
(N ~ 4 2 ) or an odd one (N 3 2 , 5 2 ). Possible systematic error of is 
argued in Section 12.51 

We turn to the analysis of the correlation-length critical exponent v. In 
Fig. |5l we plot the approximate critical exponent (pluses) 

v a (N N) = ln(Li/L 2 ) . > 

1 u 2) HdjU^/djU^N,)}]^^ 1 1 

for [2/(iVi + A" 2 )] 2 with a = V and 10 < Nt < N 2 < 20 (L 1)2 = ^N^i)- 
The parameters are the same as those of Fig. [3j Again, there emerges a 
wavy character intrinsic to the screw-boundary condition; a notable bump 
at 4 2 would be an artifact, preventing us to access to the thermodynamic 
limit systematically. The least-squares fit to these data yields v v = 1.001(93) 
in the thermodynamic limit. Uncertainty of this result is argued in Section 
1231 

2. 3. Critical behavior of the Hematic- ( quadrupolar- ) order parameter 

In this section, we investigate the critical behavior of the nematic- (quadrupolar- 
) order parameter. 
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In Fig. El we present the Binder parameter 



_ ( m Q) 

3«> 2 ' 



(12) 




nematic-order parameter is 



m 2 Q = QiQ 



(13) 



i=N/2 



The parameters are the same as those of Fig. [3j The nematic order appears 
to develop in the large- J regime, offering a sharp contrast to the VBS order 
(Fig. E]). We amended the definition of rng, Eq. f[T3"j) . so as to attain an 
intersection point of the Uq curves in Fig. otherwise, the intersection point 
disappears. Namely, we discarded the short- distance contributions within the 
rriq correlator in order to get rid of corrections to scaling. As a byproduct, 
the location of intersection point drifts significantly with respect to N; a 
subtlety of Q-based result is reconsidered in Section 12.51 

In Fig. HI we plot the approximate transition point J®(N\, N 2 ) (crosses), 
Eq. (USD, for [2/(JVi + N 2 )} 2 with 10 < N x < N 2 < 20. The parameters are 
the same as those of Fig. [3j The finite-size errors seem to be larger than 
those of . The least-squares fit to the data of Fig. H] yields an estimate 
J® = 0.151(11) in the thermodynamic limit; the error margin is considered 
in Section I2T51 

We turn to the analysis of the correlation-length critical exponent v®(Ni, N 2 ), 
Eq. (fTT|) . In Fig. [5j we plot the approximate critical exponent v®(Ni, N 2 ) 
(crosses) for [2/(JVi + A^ 2 )] 2 with 10 < Ni < N 2 < 20. The parameters 
are the same as those of Fig. [3j The least-squares fit to these data yields 

= 0.902(45) in the thermodynamic limit; a possible systematic error is 
appreciated in Section 12.51 

2.4- Fisher's critical exponent rjy,Q 

At the critical point J = J c , the quadratic moment myn obeys the power 
law, m V Q ~ 1/L 1+Vv ^, with Fisher's exponent (anomalous dimension) t]v,q 
and system size L. 

First, we consider the case of m v . In Fig. [TJ we plot the approximate 
critical exponent 



Vv(Ni,N 2 ) 



lnf 

ln(Li/L 2 ) 



-1, 



(14) 
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for [2/(JVi + N 2 )} 2 , a = V (plusses), and a = Q ( crosses) with 10 < iV\ < 
N 2 < 20. The parameters are the same as those of Fig. [3J The least- 
squares fit to these data yields r\y = 0.786(38) and i]y = —0.12(16) in the 
thermodynamic limit. 

Second, we consider Fisher's exponent for the nematic order. In Fig. [8j 
we plot the approximate critical exponent 



for [2/(JVi + iV 2 )] 2 , a — V (plusses), and a = Q (crosses) with 10 < N x < 
N 2 < 20. The parameters are the same as those of Fig. [3j The least- 
squares fit to these data yields rjq = —0.141(94) and t]q = —0.089(18) in the 
thermodynamic limit. 

2.5. Extrapolation errors of J c , v and t}v,q 

In this section, we consider possible systematic errors of J c , v, rjv, and 
t}q obtained in Figs. EU EJ and [HJ respectively. 

First, we consider the critical point J c (— J^' Q )- In Fig. HI we made 
independent extrapolations, which yield = 0.236(18) and J® = 0.151(11). 
The discrepancy (systematic error) between these results, ~ 0.09, is larger 
than the insystematic (statistical) errors, O(10 -2 ). Hence, we consider the 
former as the main source of the extrapolation error. Taking a mean value 
of and J®, we estimate J c = 0.19(9). We address a number of remarks. 
First, this estimate is consistent with the preceding one J r . — 0.285(5) [l9[ 



within the error margin; note that the preceding estimate [19j was calculated 
through the scaling of the excitation gap. Second, the result would be 
more reliable than J^; actually, the slope of J^(Ni, N 2 ) in Fig. @] appears to 
be smaller than that of J®, suggesting that the former extrapolation would 
be trustworthy. Last, the extrapolated critical point J c does not affect the 
subsequent analyses; rather, the approximate critical point J^(N\, N2) was 
fed into the formulas for the critical indices, Eqs. (ITT]). f[T4|) . and f lT5|) . 

Second, we turn to the correlation-length critical exponent v{= i> v >Q). 
In Fig. O two independent extrapolations yield v v = 1.001(93) and = 
0.902(45). The discrepancy, ~ 0.1, is comparable with the statistical error, 
~ 0.1. Taking a mean value, we obtain an estimate 



Vq(Ni,N 2 ) 



\n[m 2 Q (N 1 )/m 2 Q (N2)]\ j=j? {Ni ,n 2 ) 
HL1/L2) 



- 1 



(15) 




v = 0.95(14). 



(16) 
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Here, the error margin comes from 0.14(« VO.l 2 + 0.1 2 ) (propagation of 
uncertainty) with systematic (~ 0.1) and insystematic (~ 0.1) errors. The 
estimate, Eq. (flEj) . agrees with the preceding one v = 0.92(10) Our 
result may support the deconfmed-criticality scenario, suggesting that the 
exponent v acquires an enhancement, as compared with that of the three- 
dimensional Heisenberg universality class, v = 0.7112(5) j32|. 

Last, we consider Fisher's exponent t]v,q- m Fig. [7J the results, r\y = 
0.786(38) and rjy = —0.12(16), are obtained. The discrepancy between them, 
~ 0.9, dominates the statistical error, ~ 0.03. Considering the former as a 
error margin, we obtain rjy = 0.3(9). On the one hand, in Fig. El the 
discrepancy between rjq = —0.141(94) and r/g = —0.089(18) is negligible. 
Considering the statistical error ~ 0.1 as a main source of uncertainty, we 
arrive at r]q = —0.1(1). 

This is a good position to address a remark. As mentioned above, the 
VBS-order-based results, and rjyn, are more reliable than the nematic- 

order-based ones, J® and Vv,q- I n g enera l ; the discrete symmetry, m v , is 
more robust than the continuous one, nig, allowing us to make a system- 
atic scaling analysis even for the system size tractable with the numeri- 
cal diagonalization method. In fact, as claimed in Section III F 2 of Ref. 



33J, the numerical diagonalization method up to N = 4 x 4 is incapable 
of providing a conclusive evidence for the spontaneous magnetization of the 
two-dimensional Heisenberg antiferromagnet. Hence, tentatively, we discard 
the nematic-order-based results, and refer to the VBS-order-based ones to 
obtain crude estimates, rjy ~ 0.8 and rjQ w —0.1. These conclusions are 



comparable with the large-scale-Monte-Carlo results 29], r\y = 0.74(10) and 
r] = —0.31(6), for the Q v model; here, we made use of the scaling relation, 
1 + n = 2(3 /v, to evaluate riv and rj from (/V,/3) = [0.68(3), 0.27(2)] and 



v = 0.78(3) [29(. Possibly [29(, the exponent rj suffers from "drift" (scal- 
ing corrections), and may restore rj > through taking into account of yet 
unidentified scaling corrections. As mentioned in the Introduction, in the 



preceeding study of the S = 1-spin model [17j, an estimate i]q = 0.25(20) 
was reported. It is suggested that the anomalous dimension for the nematic 
order i]q would be almost vanishing. 

3. Summary and discussions 

The phase transition separating the VBS and nematic phases (Fig. [1]), the 
so-called deconfined criticality, was investigated for the 5* = 1-spin model on 
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the spatially anisotropic triangular lattice with the biquadratic interaction, 
Eq. (TO . So far ; the criticality has been analyzed via the scaling of the first 



excitation gap [191 ] . In this paper, evaluating both VBS- and nematic-order 
parameters, we made complementary approaches to this criticality. As a 
result, we estimate the correlation-length critical exponent v = 0.95(14) 
u v ' Q ); this estimate agrees with v = 0.92(10) [l9j], supporting the deconfined- 
criticality scenario. 

Encouraged by this finding, we put forward the analysis of Fisher's ex- 
ponent T]yQ for each order parameter. As overviewed in the Introduction, so 
far, the exponent t]q = 0.25(20) has been reported as to the S = 1-spin 
model; corrections to scaling did not admit to the estimation of r/y. We 
obtained a set of indices (?7i/,^q) = [0.3(9), —0.1(1)]; again, the exponent 
7]v suffers from scaling corrections. Alternatively, setting a = V tentatively 
(Sec. I2.5p . we arrive at crude results (^v^q) ~ (0.8,-0.1). These results 



might be reminiscent of the large-scale-simulation results 29j, v = 0.78(3), 



7] V = 0.74(10), and r\ = —0.31(6), for the Q v model. We suspect that the ex- 
ponent t)q would almost vanish through taking into account of corrections to 
scaling properly. In Ref. 34|, Fisher's exponent was estimated rather accu- 



rately through scrutinizing the local-moment distribution around a magnetic 
impurity; this idea has a potential applicability to a wide class of systems. 
This problem will be addressed in future study. 
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Figure 1: A schematic phase diagram for the S = 1-spin model on the spatially anisotropic 
triangular lattice with the biquadratic interaction, E g. (HI ), is presented. The limiting cases 
J = and J — > oo were studied in Refs. [20[ and |2lf/respectivelv. The solid (dashed) 
lines stand for the phase boundaries of discontinuous (continuous) character. The dotted 
lines are ambiguous. We investigate the phase boundary separating the nematic and VBS 
phases. 
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Figure 2: (a) We consider a spatially anisotropic triangular lattice; the Hamiltonian is 
given by Eq. ([T]). The interaction J interpolate the one- and two-dimensional lattice 
structures in the limiting cases of J — and J — > oo, respectively, (b) In order to take 
into account such a geometrical character, we implement the screw-boundary condition. 
As shown in the drawing, a basic structure of the cluster is an alignment of spins {Sj} 
(i < N). Thereby, the dimensionality is lifted to d — 2 by the bridges over the (viV)-th 
neighbor pairs through the J bonds. Technical details are explained in Section I2TT1 




Figure 3: The Binder parameter for the VBS-order parameter Uv © is plotted for 
various J and N = (+) 10, (x) 12, (*) 14, (□) 16, (■) 18, and (o) 20. The parameter j 
is set to j — 0.5. (J' is the unit of energy.) The VBS order develops in the small- J side. 
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Figure 4: The approximate critical point J?(Ni,N 2 ) ([TDJ) is plotted for [2/(iVi + N 2 )} 2 
(10 < Ni < N2 < 20) with a — V (plusses) and a = Q (crosses). The parameters are the 
same as those of Fig. [3] The least-squares fit to these data yields JY — 0.236(18) and 
J® = 0.151(11), respectively, in the thermodynamic limit. 
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Figure 5: The approximate critical exponent v a (Ni 1 N2) (fTTj) is plotted for [2/ (N1+N2)] 2 
(10 < Ni < N2 < 20) with a = V (plusses) and a = Q (crosses). The parameters are 
the same as those of Fig. [31 The least-squares fit to these data yields v = 1.001(93) and 
i>Q = 0.902(45), respectively, in the thermodynamic limit. 
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Figure 6: The Binder parameter for the nematic- (quadrupolar-) order parameter Uq 
(23) is plotted for various J and N = (+) 10, (x) 12, (*) 14, (□) 16, (■) 18, and (o) 
20. The parameters are the same as those of Fig. [3J The nematic order develops in the 
large- J side. 
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Figure 7: The approximate critical exponent T]%-(Ni,N 2 ) (QJ) is plotted for [2/(7Vi + A^ 2 )] 2 
(10 < Ni < N2 < 20) with a = V (pluses) and a = Q (crosses). The parameters are the 
same as those of Fig. [3l The least-squares fit to these data yields i]y — 0.786(38) and 
rjy = —0.12(16), respectively, in the thermodynamic limit. 
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Figure 8: The approximate critical exponent t)q(Ni, N 2 ) (|15l) is plotted for [2/(Ni + N 2 )] 2 
(10 < Ni < N2 < 20) with a = V (pluses) and a = Q (crosses). The parameters are the 
same as those of Fig. [3] The least-squares fit to these data yields t]q = —0.141(94) and 

t]q = —0.089(18), respectively, in the thermodynamic limit. 

[8] V. N. Kotov, D.-X. Yao, A. H. Castro Neto, and D. K. Campbell, Phys. 
Rev. B 82 (2010) 136402. 

[9] L. Isaev, G. Ortiz, and J. Dukelsky, J. Phys.: Condens. Matter 22 (2010) 
016006. 

[10] A. Kuklov, N. Prokof'ev, and B. Svistunov, Phys. Rev. Lett. 93 (2004) 
230402. 

[11] A.B. Kuklov, M. Matsumoto, N.V. Prokof'ev, B.V. Svistunov, and M. 
Troyer, Phys. Rev. Lett. 101 (2008) 050405. 

[12] F.-J. Jiang, M. Nyfeler, S. Chandrasekharan, and U.-J. Wiese, J. Stat. 
Mech. (2008) P02009. 

[13] K. Kriiger and S. Scheidl, Europhys. Lett. 74 (2006) 896. 

[14] J. Sirker, Z. Weihong, O.P. Sushkov, and J. Oitmaa, Phys. Rev. B 73 
(2006) 184420. 

[15] D. Poilblanc, A. Lauchli, M. Mambrini, and F. Mila, Phys. Rev. B 73 
(2006) 100403(R). 



15 



[16] Y. Nishiyama, Phys. Rev. B 85 (2012) 014403. 

[17] K. Harada, N. Kawashima, and M. Troyer, J. Phys. Soc. Japan 76 (2007) 
013703. 

[18] Y. Nishiyama, Phys. Rev. B 79 (2009) 054425. 

[19] Y. Nishiyama, Phys. Rev. B 83 (2011) 054417. 

[20] G. Fath and J. Solyom, Phys. Rev. B 51 (1995) 3620. 

[21] K. Harada and N. Kawashima, Phys. Rev. B 65 (2002) 052403. 

[22] M. A. Novotny, J. Appl. Phys. 67 (1992) 5448. 

[23] Y. Nishiyama, Phys. Rev. E 78 (2008) 021135. 

[24] N. Read and S. Sachdev, Phys. Rev. Lett. 62 (1989) 1694. 

[25] N. Kawashima and Y. Tanabe, Phys. Rev. Lett. 98 (2007) 057202. 

[26] J. Lou, A. W. Sandvik, and N. Kawashima, Phys. Rev. B 80 (2009) 
180414(R). 

[27] R. K. Kaul, Phys. Rev. B 84 (2011) 054407. 

[28] K. S. D. Beach, F. Alet, M. Mambrini, and S. Capponi, Phys. Rev. B 
80 (2009) 184401. 

[29] J. Lou and A. W. Sandvik, Phys. Rev. B 80 (2009) 212406. 

[30] T. Grover and T. Senthil, Phys. Rev. Lett. 107 (2011) 077203. 

[31] O. I. Motrunich and A. Vishwanath, Phys. Rev. B 70 (2004) 075104. 

[32] M. Campostrini, M. Hasenbusch, A. Pelissetto, P. Rossi, and E. Vicari, 
Phys. Rev. B 65 (2002) 144520. 

[33] E. Manousakis, Rev. Mod. Phys. 63 (1991) I. 

[34] A. Banerjee, K. Damle, and F. Alet, Phys. Rev. B 82 (2010) 155139. 



16 



